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A decay width calculation for a hybrid exotic meson h, with J^*^ = 1 '", is presented for the 
channel h — > ttoi. This quenched lattice QCD simulation employs Liischer's finite box method. 
Operators coupling to the h and yrai states are used at various levels of smearing and fuzzing, and 
at four quark masses. Eigenvalues of the corresponding correlation matrices yield energy spectra 
that determine scattering phase shifts for a discrete set of relative yrai momenta. Although the 
phase shift data is sparse, fits to a Breit-Wigner model are attempted, resulting in a decay width 
of about 60 MeV when averaged over two lattice sizes. 



o 
o 

(N 

in 

> 
in 
o 
o 
^ 
o 

o 

I 



X 



PACS numbers: 12.38.Gc, 13.25.-k 



I. INTRODUCTION 

Hybrid mesons are quark-antiquark pairs having va- 
lence gluons as a structural component. In some cases 
their quantum numbers are not accessible with quark 
models, and they are therefore called exotic. Examples of 
these exotics are the J^'~^ ~ 0"' , 1 ^, and 2^ mesons. 
Because these mesons contain valence gluons their verifi- 
cation is one of the signature tests of Quantum Chromo- 
dynamics (QCD). 

Efforts to determine properties of these hybrid exotic 
states are unsettled from both experimental and theoret- 
ical viewpoints 0, BS- The experimental efforts date 
back over a decade H, 0, IE S H S ^| j and currently, 
considerable resources are being devoted to their future 
study. The Jefferson Lab GlueX experimental program is 
designed to investigate exotic states. Hybrid meson stud- 
ies are also part of the COMPASS experiment at CERN 
and the CLEOc program at Cornell. 

Attempts to calculate decay widths of hybrid mesons 
have been made using the bag model the quark 

model extended by gluon flux tube degrees of freedom 
[11 m and in lattice QCD [Hill. 

Calculating properties of resonances using Euclidean 
lattice QCD simulations is not straightforward for a 
variety of reasons. On a finite lattice all states are 
bound, and furthermore, the lattice total energy of a 
two-hadron state in a decay channel is typically larger 
than the energy of the original hadron thus preventing 
decay. Aspects of these points have been discussed by 
Michael [ll[l3, DeGrand HI, Liischer H^, and by Lel- 
louch and Liischer poj . 

Thus, lattice work on hybrid mesons has concen- 
trated mainly on their mass spectrum. It can be 
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roughly classified in terms of heavy quark systems using 
static quarks with (excited) glue treated in the Born- 
Oppenheimer approximation |2ll |23| , using non- 
relativistic QCD 113, and studies using ac- 

tions with both quenched and unquenched quark dynam- 
ics 27, 28, 29, 30, 31, 32, 33, i^. Lattice work on hadron 
resonances, of which there is very little at this time, has 
recently been reviewed by Michael [T7l |. 

Hybrid meson decay widths have been studied on the 
lattice in the heavy quark limit and, recently, for light 
quarks TT] . In those studies the time dependence (slope) 
of a normalized transition matrix element computed on 
the lattice is related to a decay width via Fermi's golden 
rule [3^1. For this approach to work the lattice parame- 
ters have to be such that the resonance mass comes out 
close to the threshhold of the decay channel. 

We here choose to extract decay widths using Liischer's 
finite box method [T^ l36L l3^ . In principle Liischer's 
method is rigorous: The two-particle energy spectrum 
in a finite periodic box is related to continuum elastic 
scattering amplitudes. The spectra allow calculation of 
the scattering phase shifts at a discrete set of momenta 
owing to exact formulae derived by Liischer. A decay 
width can then be extracted by fitting a Breit-Wigner 
function, provided a resonant state is actually present. 
The applicability of this method to extract scattering 
phase shifts has been demonstrated for the 0(3) non- 
linear sigma model in 1-1- 1 dimensions l38[. the 0(4) non- 
linear sigma model in 3-1-1 dimensions j39l |4C| , nieson- 
meson scattering in 2-1-1 dimensions using QED 14 ll . and 
resonance scattering of two coupled Ising systems p2 . l43| . 

AppHcation of Liischer's method to our desired goal re- 
quires a set of operators that couple to the hybrid meson 
state and to appropriate two-meson systems matching a 
decay channel. The exotic 1 ^ meson can decay into tt&i, 
tt/i, and Trai, which are relative S-wave channels. Other 
decays are possible, but those involve relative P-waves, 
where the relative momentum is at least 27r/L, and thus 
give rise to large upward energy shifts which makes the 
simulation more difficult 18, 35, 44J. We note that ex- 
perimentally [4^ the three mesons 61 (1235), ai(1260), 
/i(1285) are close in mass, and that the oi is a vector 
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meson with C = + so that it combines naturally with 
C = of the (neutral) pion to the required charge con- 
jugation, C = +, of the exotic meson. Therefore, in this 
study, we model our two-hadron operator after the ttoi 
decay channel. 



II. LATTICE PARAMETERS 

Obtaining excited state spectra from correlation matri- 
ces that involve two-hadron operators presents a numeri- 
cal challenge in itself. On top of this more analysis steps 
are required to obtain a decay width, which essentially 
amounts to extracting a "derivative" quantity from the 
simulation. While statistical errors, and also discretiza- 
tion errors to some extent, can be easily controlled this is 
not true for the systematic uncertainties indigenous to a 
task of this kind. For this reason we have refrained from 
using large computing resources and therefore employed 
a very simple lattice action and moderate lattice sizes 
and pion masses. 

Simulations were performed using the Wilson gauge 
field action and Wilson fermions in quenched approxi- 
mation on anisotropic lattices. We will present results 
from 200 gauge configurations on 12'^ x 24 and 10'^ x 24 
lattices. 

The definition of the coupling parameters used for the 
Wilson gauge action 



E E 

X l<p<iy<4 



l-iRcTr(C/^,(x))) (1) 



is given by 



(2) 



where fiv denotes plaquette planes and oi = 02 = 
03 =: Os is the spatial and 04 =: at the temporal lat- 
tice constant. We have chosen the (bare) anisotropy 
^ = Os/ot = 2 and /? = 6.15 for the global coupling 
parameter. Parameters for the anisotropic fermion ma- 
trix 



(3) 



Q{x, y) = l S^.y - ^ f (1 - 7^)t/^(x)4.+^,y 

+ (l+7M)C^^(2/)'5r.,a+A 

are given by 



(4) 



where k is a global hopping parameter. With these con- 
ventions the relation of the latter to the (bare) Wilson 
quark mass parameter is 



A ax 



8K + Ea^)' 



(5) 



which identifies Kc = 0.125 as the critical value. 

Correlation functions for mesons from standard local 
operators tt ^ 75, p ^ 7^, and oi 7i75, i — 1,2,3 
were constructed as a matter of course employing three 
iterations of quark field smearing and gauge field 
fuzzing Also, within this setting, we adopted the 

hybrid meson operator proposed in '48l with magnetic 
type gluons 



Oh{t) 



E E 

l<i<j<3 X 



da{xtHub{xt) [i^f (ft) 



(6) 

where o, b denote color indices and Fij (x) is a product of 
SU(3) link matrices arranged in a clover pattern 

F^.ix) = U^{x)U,{x + fi)Ul{x + i^Plix) 

+ U^{x)Ul{x - fl + £')Ul{x - fi)Uf,{x - fl) 

+ Ul{x - fi)Ul{x - /t - i')U^,{x - /t - v)Uv{x - v) 

+ Ul{x~v)U^,{x~v)U,{x + fL-v)Ul{x), (7) 

which is used in the spatial planes only. Under parity we 
have VOh{t)V~^ = —Oh{t), while for the charge neutral 
{du — > uu, dd ) version O^o (t) of lO under charge conju- 
gation the derivation of COfio{t)C~^ = —Oho{t) relies on 
the presence of F - F^, specifically CFij ix)C^^ — F*j (x). 
We have enforced this relation in the simulation: Observ- 
ing that Sg[U] = Sg[U*] the configurations [U] and [U*] 
are equally probable. Thus with each [U] in the ensemble 
of 200 configurations we also include [U*] and compute 
fermion propagators for both of those. This strategy dou- 
bles the number of fermion propagators that need to be 
computed, however, charge conjugation is now numeri- 
cally exact, and this also appears to be the reason for an 
observed noise reduction of simulation signals. 

Meson masses were obtained at four values of the hop- 
ping parameter k, see Tabled A multiple mass inverter 
[49| was used to compute propagators. The resulting 
ground state masses, coming from three smearing itera- 
tions, sources set at t = 3, and effective mass function 
fits in the range i = 6 ... 11, are also listed in Tab. ^ 



TABLE I: List of hopping parameters k and the resulting tt, 
p, ai and h meson masses in units of the temporal lattice 
constant at for lattices 12^ x 24 (upper table) and 10^ x 24 
(lower table). 



K 


atrriTr 


atnip 


atrriai 


atrrih 


0.140 


0.53(4) 


0.55(3) 


0.65(4) 


0.63(25) 


0.136 


0.64(3) 


0.65(3) 


0.75(3) 


0.75(21) 


0.132 


0.75(3) 


0.75(3) 


0.86(3) 


0.85(26) 


0.128 


0.85(4) 


0.85(3) 


0.96(3) 


0.95(24) 


0.140 


0.54(5) 


0.54(3) 


0.64(4) 


0.62(24) 


0.136 


0.65(3) 


0.62(3) 


0.75(3) 


0.71(23) 


0.132 


0.74(3) 


0.73(3) 


0.85(4) 


0.81(25) 


0.128 


0.85(4) 


0.83(3) 


0.96(3) 


0.91(22) 



In order to allow extrapolations to the small pion mass 
region it is useful to study the dependence of the com- 
puted p, oi and h masses on = atmriT^. Predictions 
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for this dependence may in principle come from chiral 
pertm'bation theory [s^l, and will depend on the baryon 
being studied. For a baryon of mass M = atm the ex- 
pression 

M ^po+P2 Ml + p3 Ml + p4 Ml + q Ml ln(Af, ) (8) 

contains a collection of terms typical for xPT inspired 
models [^H^. In the case of the hybrid exotic meson, 
for example, the authors of fs^ retain only the even poly- 
nomial in JHl- For the ttoi decay channel, which is mostly 
relevant in this work, no predictions for the dependence 
of the spectral masses, say W, on M^r are available. A 
three parameter model that reflects features of JHJ is 



W = p + qx + r ln(l -|- x) with x — {atniT^y 



(9) 



The logarithmic term is purely heuristic. Its role is to 
provide curvature to the model, just like the last three 
terms in ^ do while vanishing as M^^ 0. As it turns 
out this model yields fits that are on average optimal on 
our spectral data for the combined h and ttoi systems. 
Replacing the logarithmic term in (j^ with x^/^ yields 
nearly identical results. 

We consistently use (jnj to fit all masses emerging from 
the simulation. Examples for the mesons listed in Tab. d 
are shown in Fig. ^ 

In particular the upper panel of Fig. ^ exhibits a level 
crossing between the hybrid meson mass and the n + ai 
mass, assuming a relative S-wave for the latter. The 
level crossing emerges near a: ~ which is only reached 
through extrapolation. The lower panel illustrates the 
effect of P-wave vs. S-wave decay on the lattice. The mass 
of the TT -(- TT system is shown with pious having lattice 
momenta ±27r/ (a^L). Clearly a level crossing with the p 
meson mass is harder to achieve. 

The extrapolated p meson mass, at a; = 0, shall be used 
to set the physical mass or length scale for this simula- 
tion. We obtain a* = 0.33(5) GeV"^ = 0.07(1) fm {nip = 
776 Me V). If the ai meson is used instead the scale is 
at = 0.30(3) GeV"^ = 0.06(1) fm {ma, = 1230 Me V). 
The above numbers are based on the 12^^ x 24 lattice. 
Unless otherwise indicated we will quote results using 
the p meson to set the scale. 

As a sideline it is interesting to note that the level 
crossing seen in Fig. ^ using both of the above scales, 
thus occurs within 1.35-1.49GeV, which overlaps with 
the experimental mass of the 7r(1400) resonance, accord- 
ing to [43. Indeed the 7r(1400) has the quantum numbers 
1 ^ of the hybrid exotic meson. This observation coin- 
cides with the findings of [ssf. 



III. CORRELATION MATRIX 

The description of /i — > tt -I- ai requires an operator 
for the two-meson decay channel with suitable quantum 




FIG. 1: Combinations of masses, as indicated, obtained from 
single meson operators versus the squared pion mass x = 
{atrriTr)^ and fits with the model Q. The relative S-wave 7v + 
ai mass (upper panel) reveals a level crossing near zero pion 
mass through extrapolation. The n + n mass is shown with 
pions having lattice momenta ±27r/(asL). The extrapolation 
of the p meson mass to a; = is used to set the physical scale. 



numbers. Consider 

0^+a»-k,r{'t) X! X! ^^-V,r da{xt)-f^Ua{xt) 
X y 

[db{yt)lblkdb{yt) + ub{yt)^^^k ub{yt)] (10) 

where = 1, 2, 3. The relative distance r may in principle 
be used to construct operators that transform according 
to an irrep of the hypercubic group. However, the sim- 
plest choice r = already leads to a viable operator. 
Summing over all spatial directions we thus adopt 



(11) 



fc=l 



for this simulation. 

The operators ((BJ and (|11() are the basis for calculating 
correlation functions 

CxY{t,to) = {Ox{t)Ol{to)) - {Ox{t)){0l.{to)) . (12) 
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Here X and Y stand for h or Trai, and thus establish a 
2x2 correlation matrix. The separable terms in 112|) 
are zero because of the quark flavor assignment in h 
and TTfli. The remaining (non-separable) terms in l|12|) 
contain contractions between quark fields at equal times 
when worked out with Wick's theorem. For example, 
showing flavor structure only, 

C/i,7rai [t-, to) ~ {{du)t{ddud + uuud)to ) (13) 

with time arguments t and as indicated. Equal-time 
contractions :dd: and :uu: occur only at the source 
time slice to. The corresponding propagator elements 
Q~^{xto,yto) are calculated by default. This is different 
for 

C'lrai.Trai (t, to) ~ {{dudd + duuu)t{ddud + uuud)tg) (14) 

where we encounter equal time contractions :dd: and :uu: 
at t > to- The computation of Q~^{xt,yt) is very re- 
source intensive and, if stochastic estimation is used [53| . 
then it contributes additional noise. 

We shall now argue that this problem can be circum- 
vented: At nid = rriu the contractions :dd: and :uu: give 
rise to the same propagator elements_Qj^(a;t, yt). Thus, 
replacing the (. . .)t term in H14|) by {2dudd)t and reinstat- 
ing the 7-matrices from (IIOII we observe that dd ~ d'y^jkd 
couples to oi and /i mesons. Their masses are close how- 
ever, 1230MeV and 1282MeV respectively Invoking 
a similar argument, altering the quark flavor d — > s in 
the above operator entails dj^jkd — > ^757^5 ~ Ki and, 
again, should not significantly alter the mass spectra be- 
cause the Ki meson mass of 1270MeV again is close 
to that of the ai meson. In terms of (|14|l the effect is 

C^ai,Trai(i,io) ^ {{2duds)t{2sdud)to) , 

which now has no equal time contractions, but other- 
wise is not different from 114|) when worked out in terms 
of quark propagators. Hence, dropping equal time con- 
tractions in H14|l should have little effect on the mass 
spectrum of the ttoi system, and ultimately, on the re- 
sulting decay width, because Liischer's method for com- 
puting scattering phase shifts exclusively relies on the 
mass spectra, in a finite box. 

We emphasize that the correlator element (|14|l is 
worked out using the quark flavor structure exactly as 
it emerges from (|10ll except that equal time contractions 
are neglected. All other matrix elements are not effected. 

Finally, we do not explicitly compute C7rai.h(t, to) but 
rather infer it from the hermiticity of the correlation ma- 
trix. 



IV. ANALYSIS 

For every operator used in this simulation up to three 
iterations of quark field smearing and gauge field 



fuzzing [43 were employed, using 2.5 as a strength pa- 
rameter in both cases. The 2x2 correlation matrix l|12|) 
thus expands to size 6x6, 



C(t, to) 



Ch{},h{}{t,to) Ch{},Trai{}it,to) 
C'ttoi {},/i{} (^: ^0) C'ttoi {},7rai {}(*> ^0) 



(15) 



where the entries are 3x3 matrices with elements 
Cx{k},Y{e}{tjto) built from operators Ox{fe}(^i^o)j etc, 
with k = 1, 2, 3 levels of fuzzing and smearing. The lat- 
ter is done identically at both source and sink, and thus 
the matrix C(t, to) is hermitian by construction. 

A standard analysis method is based on solving the 
generalized eigenvalue problem [s^ 



C(t,to)*(t) =C(ti,to)«'(t)A(t) 



(16) 



where ti is fixed, ^'(t) is an x matrix, its columns 
being the generalized eigenvectors, and A(t) is real diago- 
nal. We further require that C(ti, to) be positive definite. 
To ensure the latter ti should be an "early" time slice, 
here we use ti — to = 4. The generalized eigenvalue prob- 
lem (|16ll can then be cast into an ordinary one by first 
diagonalizing 

C(ti,to) = V{ti,to)D{h,to)vHh,to), (17) 

where V^(ti,to) is unitary and I?(ti,to) is real diagonal 
and positive definite. Inserting lfT7|) into ifTCjl we are lead 
to define 



C(t) = 



1 



Ft(ti,to)C(t,to)F(ti,to) 



1 



(18) 

It's eigenvalues are the same as those of the generalized 
problem lfTB|l . A(t) — diag(Ai(t) . . . AN(t)). Constructing 
C(t) merely amounts to a linear transformation among 
the set of operators Ox{fc}(ij*o) that define the correla- 
tion matrix. Because of C'(ti) = 1, the transformed set 
of operators create quantum states that are orthogonal 
and normalized at t = ti . Provided that these match the 
"true" states of the theory, the eigenvectors of C'(t) will 
stay orthogonal as t > ti increases to the extent allowed 
by the errors of the simulation. Consequently C(t) will 
be diagonal dominated for t > ti. 

From an analysis point of view there are now two op- 
tions: One may diagonalize C(t) on each time slice t sepa- 
rately and thus obtain the eigenvalues A„(t), n = 1 . . . A^, 
as proposed in [3^] . Alternatively, the diagonal elements 
of C(t) may be taken as an approximation to its eigen- 
values. 



An(t) ~ Cnn{t) 



(19) 



which involves projecting C(t) into the eigenspaces at 
fixed time slice ti, see (I18|) . The latter approach has 
the advantage that statistical fluctuations are reduced, 
see Fig. |21 for a comparison. This is plausible because 
fluctuations of the eigenvector components are effectively 



5 





2 









_9 
z 




-4 




-6 








-8 








-10 




-12 




14 




-16 




K = .1400 



10^x24 



2^0 




FIG. 2: Comparison of eigenvalues obtained from diagonaliz- 
ing C{t) on every time slice (upper panel) versus projecting 
to time slice t = 5 (lower panel), with io = 1- Results are for 
the 10^ X 24 lattice at the lightest pion mass. Fluctuations 
are much reduced using the projection technique. 



frozen. It also has the advantage to tag the eigenvalues 
to a specific eigenvector, which is important for tracking 
the quark mass dependence of the spectral levels. 

Effective mass function plateaus typically develop in 
the time interval 5 < i — to ^ 10, or so. In this region the 
projection technique yields more stable results, particu- 
larly for the excited states. An example for the 12^^ x 24 
lattice at K = 0.140 is shown in Fig.O 

Thus we will continue our analysis with the projected 
correlators and, for simplicity, refer to Cnn(t) as eigen- 
values Xn{t). Those then give rise to the spectral energies 



1 . . .6, listed in Tah.m 



Fits to those spectra with the model Q are shown in 
Fig. 0J This figure sheds light on the volume dependence 
of the spectral levels. Evidently the ground state mass 
is relatively stable against changing the lattice volume. 
On the other hand the effect on excited states is clearly 
significant, even to the extent that level crossing patterns 
differ for some of the states. This should not be surpris- 
ing because excited levels are likely to describe two-meson 
states, which are spatially large. Nevertheless, anticipat- 



- K = .1400 



1.21(9) 



0.5 



0.0 



.69(11) 



12'x24 



4 5 6 7 8 9 10 11 12 
t 

FIG. 3: Example of effective mass functions for two (pro- 
jected) eigenvalues of the correlation matrix C{t) on the 
12^ X 24 lattice at the lightest pion mass. The higher mass 
comes from X^it) and the lower mass from X2{t). 



TABLE II: Energy spectra W„, n = 1 . . . 6 from the 12^ x 24 
lattice (upper table) and the 10'^ x 24 lattice (lower table) at 
four pion masses. 



K 


0.140 


0.136 


0.132 


0.128 


atWi 


0.55(14) 


0.69(2) 


0.74(4) 


0.82(10) 


atW2 


0.69(11) 


0.79(3) 


0.89(4) 


0.93(9) 


atW-A 


0.81(10) 


0.81(4) 


0.91(4) 


1.03(7) 


atWA 


0.88(13) 


1.03(11) 


1.18(10) 


1.26(5) 


atWz 


1.21(9) 


1.13(3) 


1.25(3) 


1.31(4) 


atWe 


0.96(5) 


1.42(4) 


1.62(3) 


1.85(4) 


atWi 


0.57(7) 


0.69(11) 


0.75(11) 


0.78(6) 


atW2 


0.59(11) 


0.71(5) 


0.78(4) 


0.81(5) 


atW-i 


0.61(6) 


0.74(5) 


0.80(4) 


0.91(5) 


atWi 


0.91(11) 


1.12(8) 


1.14(13) 


1.24(9) 


atWs 


1.20(9) 


1.42(11) 


1.64(9) 


1.78(9) 


atWe 


1.31(6) 


1.43(6) 


1.55(9) 


1.62(13) 



ing results, the volume effect on the scattering phase shift 
ultimately turns out to be only moderate. 

Another comment on Fig. ^ is that, although the 6x6 
correlation matrix gives rise to six eigenvalues, the num- 
ber of physical states on the lattice is likely to be a lesser 
number because, typically, hadronic level spacings are 
of the order of a few hundred MeV. Thus we entertain 
the possibility that the group of the three lower levels in 
Fig. ^describe the same state, the ground state, whereas 
the upper levels belong to two-meson states with some 
degree of interaction energy due to their relative motion. 
This point will become more plausible in terms of the 
corresponding scattering phase shifts. 
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0.5 



10"x24 



0.0 ' 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
(am r 

FIG. 4: Mass spectra obtained from the 12^ x 24 and 10^ x 24 
lattices versus x = (atm^)^. Error bars are omitted for clarity 
except for one level where the errors shown are typical of all 
levels. Fits with the model @ are shown as lines. 



V. SCATTERING PHASE SHIFTS 

The spectral energies Wn computed on the lattice give 
rise to a discrete set of relative tt+oi momenta fc„ by 
solving the relativistic dispersion relation 



(20) 



Note that the resulting momenta fc„ relate to spectral 
masses and thus are continuous numbers (not subject to 
lattice discretization). Only those levels Wn will be used 
that fall within the elastic region, 



(m^ + nia^ ) <Wn <2 {rriT, + ) . 



(21) 



Continuum S-wave scattering phase shifts Sa{k) are then 
computed at a discrete set of momenta /c„ using Liischer's 
formula 



tan(5o(fc„ 



277 



(22) 



Here Z{l\q^) is a generalized C-function, and Lg = Lag 
is the physical size of the spatial box, using the bare 
anisotropy as = 2at- 

If the number of available data points is sufficient , then 
one may attempt a fit to a Breit-Wigner function [Saj . 



tan 5i){k) 
where W{k) 



r/2 



So - W{k) 



fc2 



fc2 



(23) 
(24) 



The resonance energy Eq and the decay width F are fit 
parameters. However, a successful fit can only be ex- 
pected if the underlying physics indeed supports an iso- 
lated resonance. Such a fit actually fails for all spectra 
computed at the four pion masses, or rather x — (atTO^)^, 
as they appear in Fig. 01 This is not be surprising be- 
cause those data points are far away from a level cross- 
ing between the h and the 7r-|-ai masses as evident from 
Fig. n It is necessary to extrapolate the spectral masses 
to a; = near the level crossing. The model © has been 
used for this purpose. We present the extrapolated spec- 
tra in Tab, mil along with the corresponding momenta A;„ 
and scattering phase shifts (5o(A;,i) for those levels which 
fall into the elastic region H21|l . 



TABLE III: Extrapolated energy spectra W„ using the model 
@, resulting momenta fe„, and S-wave scattering phase shifts 
on lattices 12^ x 24 (upper table) and 10^ x 24 (lower table). 
Missing entries for atfe„ and (5o(fen) correspond to energy levels 
outside of the elastic region 12H . 



n 


atWn 


atkn 


5o(fc„) 


1 


0.64(16) 


0.16(1) 


74.6(11.6) 


2 


0.63(25) 


0.15(2) 


65.8(14.1) 


3 


0.19(5) 






4 


0.47(11) 


0.03(1) 


2.1(1.8) 


5 


0.31(11) 






6 


0.26(5) 






1 


0.66(21) 


0.20(2) 


77.1(10.5) 


2 


0.61(3) 


0.16(1) 


50.0(6.1) 


3 


0.49(4) 


0.06(1) 


6.0(3.1) 


4 


0.44(10) 


0.01(1) 


0.2(0.5) 


5 


0.36(10) 






6 


0.33(4) 







The phase shift data are very sparse and do not alone 
resolve the functional form of the fit model, such as Breit- 
Wigner. In fact, attempts of Levenberg-Marquardt fits 
using (|23|l only returned stable results for the resonance 
energy parameter atEo, while the width parameter OfF, 
being an indicator for a derivative, was left undeter- 
mined due to large standard errors. Nevertheless, it is 
evident from Tab. Illll that the phase shift data are clus- 
tered around two regions of atkn, namely w 0.15-0.16 
and w 0.03 for the L = 12 lattice, and « 0.16-0.20 
and « 0.01-0.06 for the L ~ 10 lattice. This suggests 
that no more than two distinct physical states are un- 
covered by the simulation. Under this assumption the 
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data may be analyzed as follows: Denoting the weighted 
(x^ = min) averages for each of the clustered momenta 
by fci,2, the corresponding energies by uji,2 = atW(ki^2), 
and Ti.2 = tan(5o(fci,2)j we obtain a set of two equations 
from (|23() . for each lattice, which are solved exactly by 



atEo 



T1T2 



T2 



T2 - Tl 



(25) 
(26) 



The resulting parameter values are listed in Tab. IIVI 
There, the uncertainties for ajF and atEo are computed 
as follows: The statistical (jackknife) errors for atfc„, as 
they appear in Tab. IIIll give rise to errors Afci^2 for the 
weighted momentum averages ki^2- Repeating the analy- 
sis procedure described above a few thousand times with 
momenta ki^2 = ^1,2 + ?AA;i_2, where ^ is a normal dis- 
tributed random deviate with variance one, then yields 
the uncertainties given in Tab. IIVI Also the ai meson 
mass, which enters the fit model via (|24|l . was subjected 
to the same randomization. The errors given in Tabic Hvl 
are the standard deviations resulting from the random- 
ization, they are reminiscent of statistical errors. Table 
1^3 also contains the physical values for the decay widths 
and the resonance energies using the p meson to set the 
mass scale. Setting the scale with the oi meson mass, 
results in widths of 39(29)MeV and 108(48)MeV for the 
12^ X 24 and 10^ x 24 lattices respectively. 

Finally, Fig. |31 shows the combined phase shift data 
from the two lattices on a common physical scale. 
The curves are Breit-Wigner interpolations as explained 
above. 



VI. ERRORS 

All errors cited in this paper are statistical, and are 
derived from a standard jackknife procedure [56j |. The 
hybrid meson operator, involving gauge link paths, ap- 
pears to be the major source of those. 

The effects of systematic errors on the results of this 
simulation are difficult to assess. In principle this can 
only be done by repeating it with various different choices 
of lattice and analysis model parameters. Probably the 
largest sources of systematic error stem from curve fitting 
and extrapolation techniques. As a check on the extrap- 
olations we have also done fits with a x^/^ term in place 
of the logarithmic term in ©. The results. Tab. IIVI did 
not change much, within statistical errors. If only linear 



TABLE IV: Results for decay widths F and resonant energies 
Eo for two lattice sizes L in units of at, and with a physical 
scale set by the p meson mass. 



L 




atEo 


r[MeV] 


Eo[GeY] 


12 


0.012(9) 


0.63(3) 


35(26) 


1.88(8) 


10 


0.032(14) 


0.62(3) 


97(43) 


1.84(9) 



120 



60 - 








0.0 



0,3 0.6 

k[GeV] 



0.9 



FIG. 5: Scattering phase shifts So{k„) from the lattices 12"^ x 
24 (filled circles) and 10^ x 24 (open circles). The solid and 
dashed curves are Breit-Wigner interpolations according to 
l|25|l and l|26|l . The filled and open box plot symbols indicate 
the respective x^-weighted averages over data points. 



terms are retained the extrapolated masses slightly shift 
upward, ultimately resulting in a slight increase of decay 
widths on the order of « lOMeV. 

Another source of systematic error comes from postu- 
lating a Breit-Wigner model. Given the sparsity of data 
points it is inconclusive that the physical phase shifts will 
indeed follow a Breit-Wigner form. In order to resolve 
this problem the simulation would have to be repeated 
at several values of the gauge coupling f3, thus mapping 
out some sort of continuous curve So^k) vs. k. The re- 
sults of this work do rely on the a priori assumption that 
the simulation data follow a Breit-Wigner model. 

On the other hand, adopting the less stringent crite- 
rion that a resonance is present if the phase shift data 
passes through 90°, the simulation results clearly indi- 
cate the presence of such. This, in itself, is a significant 
outcome of this project. Although this does not help 
putting bounds on the systematic error of F, the results 
for the resonance energy, Eq ^ 1.9 GeV, are remarkably 
stable. A decay width, on the other hand, essentially 
comes from derivative data and as such is prone to a 
significantly larger error. 

Systematic errors are also caused by finite size effects. 
At first sight, judging by the small difference of the p me- 
son masses on the L = 12 and L = 10 lattices, see Tabled 
those appear to be small. Finite size effects should be ex- 
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pected to be much larger for larger-sized hadrons like the 
Oi for example. This is particularly true for two-hadron 
systems studied in this work. For example, the spectra 
displayed in Fig.^are significantly different, particularly 
for excited states on the L — 12 and L — 10 lattices, their 
size though being quite similar. Again these effects can 
only be studied by repeating this simulation with several 
lattices of different sizes. 

Obtaining a single scattering phase data point requires 
evaluating up to three effective mass functions - one for 
the correlator matrix eigenvalue, one for the p meson to 
set the physical scale, and one for the n meson. The 
variability in choosing which time slices of the corre- 
lation functions to use in fitting effective masses pro- 
duces a variability in the decay width. Here, one usu- 
ally wants to maximize the plateau width of the effective 
mass functions to optimize the statistical error. Reduc- 
ing the plateau width to estimate a systematic error is of 
limited value. 



VII. CONCLUSION 

Decay widths for the hybrid exotic meson with J^^' = 
f ^, calculated using Liischer's method, are in the range 
35 to 97 MeV with statistical errors of about 30 MeV 
using the p meson to set the scale. The lower value for 
the width came from using extrapolated energy spectra 
on a 12"^ X 24 lattice and the higher value came from 



using extrapolated spectra on a 10^ x 24 lattice. If the 
ai meson sets the scale, then the widths for these two 
lattices range from 39 to 108 MeV with statistical errors 
of about 40 MeV. 

The number of data points available to fit Breit-Wigner 
functions is very sparse, the reason being that many en- 
ergy levels fell outside the elastic region where phase 
shifts using Liischer's formulae cannot be computed. 
Overcoming this limitation requires use of a larger cor- 
relation matrix. This can be accomplished by adding 
more smearing levels, and spatially extended operators 
in the individual correlators, or possibly by using a cou- 
pled channel type analysis in which more than one decay 
channel is represented in the matrix. Several values of 
the coupling parameter f3 should also be employed to 
generate more phase shift data points. 

Using the p meson to set the scale, the resonance mass 
of the hybrid meson in this simulation was 1.9(1) GeV, 
and in contrast to the decay width, the resonance mass 
was well determined by the simulation. This unexpected 
result leads to a final comment that, historically, hadron 
mass calculations within lattice QCD have been done us- 
ing single-hadron operators, ignoring the fact that most 
hadrons are resonances and thus are unstable [l3|. We 
have taken this decay aspect seriously. Although the nu- 
merical values for the decay widths serve as a guide only, 
the approach of extracting hadron masses as resonance 
energies, using Liischer's method, should also be given 
serious consideration. 



[1] T. Barnes (2003), nucI-th/0303032. 
[2] T. Barnes (2000), nucl-th/0009011. 

[3] T. Barnes, Acta Phys. Polon. B31, 2545 (2000), hep- 
ph/0007296. 

[4] D. Aide et al. (IHEP-Brussels-Los Alamos- Annecy 

(LAPP)), Phys. Lett. B205, 397 (1988). 
[5] D. R. Thompson et al. (E852), Phys. Rev. Lett. 79, 1630 

(1997) , hep-ex/9705011. 

[6] A. Abele et al. (Crystal Barrel), Phys. Lett. B423, 175 

(1998) . 

[7] G. S. Adams et al. (E852), Phys. Rev. Lett. 81, 5760 

(1998) . 

[8] S. U. Chung et al. (E852), Phys. Rev. D60, 092001 

(1999) , hep-ex/9902003. 

[9] E. I. Ivanov et al. (E852), Phys. Rev. Lett. 86, 3977 

(2001), hep-ex/0101058. 
[10] A. R. Dzierba et al. (2005), hep-ex/0510068. 
[11] T. Barnes, F. E. Close, F. de Viron, and J. Weyers, Nucl. 

Phys. B224, 241 (1983). 
[12] T. Barnes, F. E. Close, and E. S. Swanson, Phys. Rev. 

D52, 5242 (1995), hep-ph/9501405. 
[13] F. E. Close and P. R. Page, Nucl. Phys. B443, 233 

(1995), hep-ph/9411301. 
[14] C. McNeile and C. Michael (UKQCD), Phys. Rev. D73, 

074506 (2006), hep-lat/0603007. 
[15] C. McNeile, C. Michael, and P. Pennanen (UKQCD), 

Phys. Rev. D65, 094505 (2002), hep-lat/0201006. 



[16] C. Michael, Nucl. Phys. B327, 515 (1989). 

[17] C. Michael, in Proceedings of The 23rd International 
Symposium on Lattice Field Theory, 25-30 July 2005, 
Trinity College, Dublin, Ireland, PoS(LAT2005)008 
(2005). 

[18] T. A. DeCrand, Phys. Rev. D43, 2296 (1991). 

[19] M. Liischer, Nucl. Phys. B364, 237 (1991). 

[20] L. Lellouch and M. Liischer, Commun. Math. Phys. 219, 

31 (2001), hep-lat/0003023. 
[21] K. J. Juge, J. Kuti, and C. J. Morningstar, Phys. Rev. 

Lett. 82, 4400 (1999), hep-ph/9902336. 
[22] K. J. Juge, J. Kuti, and C. J. Morningstar, Nucl. Phys. 

Proc. Suppl. 63, 326 (1998), hep-lat/9709131. 
[23] S. Perantonis and C. Michael, Nucl. Phys. B347, 854 

(1990). 

[24] T. Burch, K. Orginos, and D. Toussaint, Phys. Rev. D64, 

074505 (2001), hep-lat/0103025. 
[25] I. T. Drummond, N. A. Coodman, R. R. Horgan, H. P. 

Shanahan, and L. C. Storoni, Phys. Lett. B478, 151 

(2000), hep-lat/9912041. 
[26] T. Manke, L T. Drummond, R. R. Horgan, and H. P. 

Shanahan (UKQCD), Nucl. Phys. Proc. Suppl. 63, 332 

(1998) , hep-lat/9709001. 

[27] C. McNeile et al., Nucl. Phys. Proc. Suppl. 73, 264 

(1999) , hep-lat/9809087. 

[28] C. W. Bernard et al. (MILC), Phys. Rev. D56, 7039 
(1997), hep-lat/9707008. 



9 



[29] C. W. Bernard et al., Nucl. Phys. Proc. Suppl. 53, 228 

(1997), hep-lat/9607031. 
[30] P. Lacock and K. Schilling (TXL), Nucl. Phys. Proc. 

Suppl. 73, 261 (1999), hep-lat/9809022. 
[31] P. Lacock, C. Michael, P. Boyle, and P. Rowland 

(UKQCD), Nucl. Phys. Proc. Suppl. 63, 203 (1998), hep- 

lat/9708013. 

[32] P. Lacock, C. Michael, P. Boyle, and P. Row- 
land (UKQCD), Phys. Lett. B401, 308 (1997), hep- 
lat/9611011. 

[33] P. Lacock, C. Michael, P. Boyle, and P. Row- 
land (UKQCD), Phys. Rev. D54, 6997 (1996), hep- 
lat/9605025. 

[34] C. Bernard et al (MILC), Nucl. Phys. B (Proc. Suppl.) 

119, 260 (2003), hep-lat/0209097. 
[35] C. McNeile and C. Michael (UKQCD), Phys. Lett. B556, 

177 (2003), hep-lat/0212020. 
[36] M. Liischer, Nucl. Phys. B354, 531 (1991). 
[37] M. S. Cook, Exotic Meson Decay Widths using Lattice 

Quantum Chromodynamics (Ph.D. Dissertation, Florida 

International University, Miami, Florida, 2006). 
[38] M. Liischer and U. Wolff, Nucl. Phys. B339, 222 (1990). 
[39] M. Gockeler, H. A. Kastrup, J. Westphalen, and 

F. Zimmermann, Nucl. Phys. B425, 413 (1994), hep- 

lat/9402011. 

[40] F. Zimmermann, J. Westphalen, M. Gockeler, and H. A. 
Kastrup, Nucl. Phys. Proc. Suppl. 30, 879 (1993), hep- 
lat/92 11029. 

[41] H. R. Fiebig, A. Dominguez, and R. M. Woloshyn, Nucl. 

Phys. B418, 649 (1994). 
[42] C.R. Gattringer and C.B. Lang, Nucl. Phys. B 391 (1993) 



463. 

[43] C.R. Gattringer , I. Hip and C.B. Lang, Nucl.Phys. B 

(Proc. Suppl.) 30 (1993) 875. 
[44] R. D. Loft and T. A. DeCrand, Phys. Rev. D39, 2692 

(1989). 

[45] S. Eidelman et al.. Physics Letters B 592, 1+ (2004), 

URL http : //pdg . Ibl . gov 
[46] C. Alexandrou, S. Giisken, F. Jegerlehner, K. Schilling, 

and R. Sommer, Nucl. Phys. B414, 815 (1994), hep- 

lat/9211042. 

[47] C. Albanese et al., Phys. Lett. B192, 163 (1987). 

[48] C. Bernard et al, Phys. Rev. D56, 7040 (1997). 

[49] A. Frommer, B. Nockel, S. Giisken, T. Lippert, and 

K. Schilling, Int. J. Mod. Phys. C6, 627 (1995), hep- 

lat/9504020. 

[50] in Chiral Dynamics: Theory and Experiment, edited by 

A. M. Bernstein and B. R. Holstein (Springer- Verlag, 

Berlin, Heidelberg, New York, 1995), vol. 452 of Lecture 

Notes in Physics. 
[51] J. Smit, Introduction to Quantum Fields on a Lattice 

(Cambridge University Press, Cambridge, UK, 2002). 
[52] S. V. Wright et al., Nucl. Phys. Proc. Suppl 109A, 50 

(2002), hep-lat/01 11053. 
[53] J. N. Hedditch et al., Phys. Rev. D72, 114507 (2005), 

hep-lat/0509106. 
[54] C. Michael and J. Peisa (UKQCD), Phys. Rev. D58, 

034506 (1998), hep-lat/9802015. 
[55] J. R. Taylor, Scattering Theory (Wiley, New York, 1972). 
[56] B. Efron, SIAM Review 21 (1979) 460. 



